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ABSTRACT 



GRB 991014 is one of the shortest gamma-ray bursts detected so far with the Wide 
Field Cameras aboard BeppoSAX, both in 7-rays and X-rays. The duration is 9.6 s in 2- 
28 keV and 3.2 s in 40 to 700 keV (as measured between the times when 5 and 95% of the 

OO . 

{SJ , burst photons have been accumulated). We refine the Interplanetary Network annulus 

of the burst, present the detection of the X-ray afterglow of GRB 991014 within this 
\ refined annulus, and discuss X-ray and 7-ray observations of the prompt and afterglow 

emission. Except for the briefness of the prompt event, no other unusual aspects were 

1-^ ■ found in the prompt and afterglow observations as compared to such measurements in 

1 ■ previous gamma-ray bursts. 
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1. Introduction 

The Wide Field Cameras (WFCs) on board the BeppoSAX satellite have, since launch in 
April 1996, detected and localized over thirty gamma-ray bursts (GRBs). The accurate (within 
a few arcminutes) and swift (within a few hours) localizations have enabled the first-time detec- 
tions of multi-wavelength afterglow emission. The optical spectra of eight such afterglows revealed 
red-shifted imprints of identifiable spectral features with z's ranging between 0.43 (GRB 990712, 
Galama et al. 1999) and 3.42 (GRB 971214, Kulkarni et al. 1998), with one exceptional redshift 
of z = 0.0085 if GRB 980425 is associated with SN1998bw (Galama et al. 1998). 

GRB 991014 triggered the Gamma-Ray Burst Monitor on BeppoSAX at 1999, Oct. 14.911508 
UT (Tassone et al. 1999). It was simultaneously detected with WFC unit 2 at an off-axis angle 
of 15° in a total field of view of 40° x 40°. The detection was of low significance and the satellite 
attitude configuration was poorly determined so that the error radius of the GRB location was 
relatively large at 6'. The location was communicated about 9 hrs after the burst. Within one 
day after the burst, a preliminary Interplanetary Network (IPN) annulus was published (Hurley 
& Feroci 1999) which intersects the WFC error region and reduces the size of the error region by 
65%. Optical and radio follow-up observations were performed but no afterglow was detected at 
these wavelengths. The sensitivity limit in the optical i?-band is 23.1 mag for observations starting 
11.2 hrs after the burst (Thorstensen et al. 1999; Uglesich, Halpern & Thorstensen 1999) and at 
radio wavelengths 125 /jJy (8.46 GHz) and 220 /iJy (1.4 GHz) for observations starting 1.68 d after 
the burst (Taylor, Frail & Kulkarni 1999). 

GRB 991014 also triggered the Burst And Transient Source Experiment (BATSE) on the 
Compton Gamma-Ray Observatory (trigger no. 7803, Giblin, Kippen & Sahi 1999). As measured in 
50-300 keV peak flux at 64 ms resolution, the burst is in the top 32% of the BATSE bursts. In fluence 
(>25 keV), it ranks in the top 58%. The T90 duration (i.e., the time interval that encompasses 90% 
of the detected photons above the background starting from when 5% is accumulated) according 
to these data is 4.67 ± 0.47 s. The hardness ratio, as determined from the ratio of the fluence in 
BATSE channels 3 (100-300 keV nominal, fluence in erg cm" 2 ) to that in channel 2 (50-100 keV), 
is H 32 = 1.36 ±0.13. 

GRB 991014 is, within the set of WFC-detected bursts, the shortest in 7-rays which was 
followed up with sensitive X-ray observations (GRB 980326 has approximately the same duration 
as GRB 991014 but was not followed up). Nevertheless, it is unlikely that GRB 991014 is a member 
of the class of short GRBs that have been recognized for many years (Norris 1984, Hurley 1992), 
that tend to have harder energy spectra than the long ones (Dezalay et al. 1992, 1996; Kouveliotou 
et al. 1993), and from which so far no afterglow has been detected in any wavelength regime. 
Giblin et al. (1999) identify GRB 991014 as belonging to the class of long bursts based on the 
BATSE-determined duration and spectral hardness. 

We present 7-ray and X-ray measurements, both of the burst event and the afterglow of 
GRB 991014. 
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2. Instrumentation and X-ray follow-up observations 

Our GRB and X-ray afterglow measurements were carried out with the following instruments 
on board BeppoSAX (Boella et al. 1997b). The Gamma-Ray Burst Monitor (GRBM; Amati et al. 
1997 and Feroci et al. 1997) comprises 4 lateral shields of the Phoswich Detector System (PDS, 
Frontera et al. 1997) and has a bandpass of 40 to 700 keV. The normal directions of two shields are 
each co-aligned with the viewing direction of a WFC unit. The GRBM has 4 basic data products per 
shield for a GRB: a time history of the 40 to 700 keV intensity with a variable time resolution of up 
to 0.48 ms, 1 s time histories in 40 to 700 and >100 keV, and a 256-channel spectrum accumulated 
each 128 s (independently phased from GRB trigger times; 240 of these channels contain scientific 
data up to 650 keV). The WFC instrument (Jager et al. 1997) consists of two identical coded 
aperture cameras each with a field of view of 40° x 40° full-width to zero response and an angular 
resolution of about 5'. The bandpass is 2 to 28 keV. The narrow field instruments (NFI) consist of 
the low-energy - 0.1 to 10 keV - and the medium-energy - 2 to 10 keV- concentrator spectrometer 
(LECS and MECS respectively, see Parmar et al. 1997 and Boella et al. 1997a respectively), the 
PDS (15 to 300 keV), and the high-pressure gas scintillation proportional counter (4 to 120 keV, 
Manzo et al. 1997). 

GRB 991014 was declared a target of opportunity (TOO) for the NFI and a follow-up obser- 
vation was performed between October 15.4 (13 h after the burst) and October 16.4 UT (we will 
refer to this observation as TOOl). The exposure time for the MECS was 38.1 ks. Two sources 
were identified in MECS data that are within the combined WFC/IPN error region (In 't Zand et 
al. 1999a). Neither of them showed clear variability. Therefore, a second observation was carried 
out (T002) between October 24.8 and 25.9 UT with a MECS exposure time of 48.2 ks. The 
X-ray afterglow was identified with the one source that had disappeared in this second observation: 
SAX J0651.0+1136 (In 't Zand et al. 1999b). The burst is at a fairly low Galactic latitude with 
(Z // ,& // )=(202?5,+5?2). 

3. Analysis 

3.1. The burst event 

In figure 1, the time profile of the burst is shown in various bandpasses. The profile is com- 
paratively simple with a single peak at all energies. Within 1 s, the peak occurs simultaneous at 
all energies. The peak profile appears to flip from below to above 40 keV (i.e., it has a fast-rise 
exponential-decay profile below 40 keV and an exponential-rise fast-decay one above that). The 
burst lasts longest at the lowest energies. In the 2-28 keV band the T90 duration is 9.6 ± 1.3 s, 
while in 40 to 700 keV this is 3.2 ± 0.7 s. We generated a spectrum from the 7-s interval when 
the WFC signal was brightest. This spectrum can be modeled by a power law function N(E)(:)E T 
phot s _1 cm _2 kcV _1 with photon index T = —0.6 ± 0.2 assuming an interstellar absorption of 
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A^h = 2.5 x 10 21 cm -2 (according to an interpolation of the HI maps by Dickey & Lockman 1990). 
The peak flux in 1-s time resolution following this spectrum is (3.2 ±0.3) x 10 -8 erg s _1 cm -2 (2-10 
keV). We estimate the 2-10 keV fluence at 1.0 x 10 -7 erg cm -2 . For 2 to 28 keV these numbers 
should be multiplied by 4.1. 

The WFC observation on this field was ongoing for 1.15 d before the occurrence of the GRB. A 
total of 39.5 ksec net exposure time was accumulated. The GRB was not detected during this time, 
with a 3<7 flux upper limit of 1.2 x 10 -10 erg s _1 cm~ 2 (2-10 keV, assuming the same spectral shape). 
Beyond 8 s after the burst onset, only 972 s of net exposure time were left in the observation. We 
could not find any signal from the burst position on time scales of 5 s, 10 s, or 972 s, and in 2-10 
or 2-28 keV. The 3a upper limit for 972 s in 2-10 keV is 7 x 10~ 10 erg s^cm" 2 . 

In 40 to 700 keV, the peak intensity is (6.4 ± 0.5) x 10 2 c s _1 in a time resolution of 1 s 
and twice as large, (1.2 ± 0.2) x 10 3 c s _1 , in a resolution of 62.4 ms. The 1 s peak flux is 
(4.2 ± 0.5) x 10 -7 erg s~ 1 cm -2 and the fluence (9 ± 1) x 10 -7 erg cm~ 2 . This brings the X-ray 
to 7-ray fluence ratio to 11% which is relatively high though not exceptional. A Fourier power 
spectrum of the high time resolution data of the GRBM shows no significant features between 3 Hz 
and 64 Hz. 

Beyond 40 keV the burst shows no spectral evolution, although the statistical quality does not 
permit a sensitive analysis. The photon index as derived from a power law fit to the two-channel 
1 s GRBM data shows no evolution with T ranging from —2.0 ± 0.7 to —1.9 ± 0.3 over the three 
brightest time bins, but the longer duration in the WFC band with respect to the GRBM band 
suggest an overall hard-to-soft evolution. The difference in T between the GRBM and WFC data 
suggests a break in the spectrum between the WFC and GRBM bands. 

The statistical quality of the 256-channel GRBM data of GRB 991014 does not permit a 
sensitive spectral analysis. 

3.2. A refined IPN annulus 

The receipt of final Ulysses ephemeris and timing data, and of the final BeppoSAX ephemeris 
data, has allowed us to reduce the width of the Ulysses/GRBM IPN annulus from 3'.3 (Hurley & 
Feroci 1999) to l'.5 (3 sigma). The final annulus, fully contained within the old one, is shown in 
figure 2. The initial error box, formed by the intersection of the preliminary IPN annulus with 
the WFC error circle, had an area of 36 square arcminutes. The final error box, formed by the 
intersection of the refined IPN annulus with the final NFI error circle, has an area of approximately 
4 square arcminutes. 
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3.3. The X-ray afterglow 

The X-ray field around GRB 991014 is somewhat complicated. Several faint X-ray sources 
can be recognized in a circular field of diameter 25' around the WFC position. Because of this 
complication, we carried out a maximum-likelihood analysis of the MECS imaging data. In com- 
parison to the standard analysis methods, the maximum-likelihood method is able to fully take into 
account the point spread function (PSF) which 1) is critical in disentangling sources that are so 
close to each other that their PSFs overlap considerably, 2) ensures that all photons detected from 
each point source are retrieved, and 3) takes into account the Poissonian nature of the data. The 
maximum-likelihood analysis was performed in two steps. In the first step a map was generated 
of detection significance, and in the second step the flux and position of all point sources were 
extracted. For details we refer to the Appendix. The significance map for TOOl data is given in 
figure 2. 

For the combined MECS TOOl and T002 imaging data, four point sources are required to 
satisfactorily describe the data. Their positions are indicated in figure 2. The elongated structure 
within the initially reported as well as refined IPN/WFC error box is explained by two point sources 
separated by We designate the upper (most northern) source as SI and the lower one as S2. 

The photon count numbers were determined in various time intervals by maximum-likelihood 
fits of the 4 point source strengths - while keeping the positions fixed - and the background level. 
The average MECS (units 2 and 3 added) photon count rate during TOOl is (3.5 ±0.4) x 10 -3 c s _1 
for SI and (2.2 ± 0.4) x 10~ 3 c s" 1 for S2. We divided TOOl in 3 time intervals to search for 
variability of the two sources, without a definite result (the maximum difference in photon count 
rate is for SI, but only at 1.4cr). In T002 the average rates are < 1.2 x 10~ 3 c s _1 for SI (3c upper 
limit) and (2.5 ±0.3) x 10~ 3 c s _1 for S2. From these latter numbers we identify SI as the afterglow 
source of GRB 991014. SI is equal to SAX J0651.0+1136. Its position is a 200 o=6 h 51 m 02.9 s , 
^2000=±ll°36 / 03 // . The statistical error radius is 40" (99% confidence) but the systematic error 
enlarges the final error radius to The combined IPN/NFI error box is depicted in figure 2. 
The ratio of the upper limit in T002 and the average flux in TOOl implies that the decay index 
is steeper than —0.4. 

A spectrum was generated for the afterglow from the LECS and MECS TOOl data by 
maximum-likelihood fits of the 4 point source strengths and background level in four logarithmically- 
sized bands between 0.2 and 2.2 keV for the LECS data and eight logarithmically-sized bands 
between 1.6 and 10 keV for the MECS data (2 units). A simple absorbed power law model with 
Galactic interstellar absorption of iVn = 2.5 x 10 21 cm -2 (following Dickey & Lockman 1990) fits the 
data well (\ 2 = 5.94 for 6 dof) with a photon index of V = —1.53 ±0.25. The average (unabsorbed) 
flux is (3.5 ± 0.5) x 10~ 13 erg s _1 cm" 2 (2-10 keV) which is a factor of 10 5 smaller than the 2-10 
keV peak flux of the gamma-ray burst itself. This implies that the energy emitted in the X-ray 
afterglow is comparable to that in the burst itself, and a decay index between the WFC and the 
MECS fluxes of —1.0 ± 0.2 which is not extraordinary quick when compared to other cases (e.g., 
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Piro 2000). The data provide a constraint on A^: it is less than 8.6 x 10 21 cm 2 . 

No signal was seen in the PDS data of TOOL The 2a upper limit in 15-100 keV is 30 times 
the flux predicted by an extrapolation of the MECS-measured spectrum. 

The NFI field seems complicated, with a concentration of sources near the GRB position. Ac- 
tually this number of sources is not a surprise given the number of background sources expected. 
Giommi, Fiore & Perri (1999) have estimated \ogN(> S)-logS distributions for 2-10 keV back- 
ground sources in BeppoSAX observations at high Galactic latitudes and find that for a flux limit 
of S = 1.0 x 10~ 13 erg s _1 cm -2 (2-10 keV) there are approximately 10 sources per square degree. 
The chance probability to find 3 or more sources in a 25'-diameter circle is 18%. The Galactic 
latitude is fairly low for GRB 991014. Therefore, one may actually expect more sources than those 
already counted at high Galactic latitude, particularly above 2 keV given Nr = 2.5 x 10 21 cm~ 2 . 
The concentration towards the GRB may be a selection effect because the sensitivity of the MECS 
increases towards the center position due to vignetting of the concentrators. 

4. Conclusion 

The constraints that could be derived for the indices of the spectrum and temporal decay are 
not unusual as compared to the X-ray afterglows of other GRBs (for a review of X-ray afterglows, 
see for example Piro 2000). Also, the spectrum of the prompt emission is not unusual. The upper 
limit on the optical afterglow is fairly faint. This is not the result of Galactic extinction which in 
this direction is a mild Ay ~ 1.4, based on HI maps by Dickey & Lockman (1990). It may be the 
result of extinction in a host galaxy. We remark that it is not exceptional that no optical afterglow 
was detected. This applies to about half of all well-localized GRBs (e.g., Kulkarni et al. 2000). 

We conclude that GRB 991014 does not exhibit behavior in the prompt spectrum or the 
afterglow that would make it distinct from other longer bursts. 

We are grateful to the staff of the BeppoSAX Scientific Operation Center, the Mission Planning 
Team and the Science Data Center for their support in obtaining and processing the data. We thank 
Tim Giblin for useful discussions. KH is grateful for Ulysses support under JPL Contract 958056, 
and to the BeppoSAX guest investigator program and NASA grant NAG5-9126 for IPN support. 
The BeppoSAX satellite is a joint Italian and Dutch program. 
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APPENDIX 

The image analysis is based on maximum-likelihood ratio (MLR) tests performed at certain 
user-defined grid points in a sky field around the object of interest, to search for point sources on 
top of a flat background. At each grid point (x s ky,y s ky) one determines the maximum likelihood 
under two hypotheses: Ho is a description of the data in terms of the flat background model only 
and Hi a description in terms of the flat background model and a point source at (x s k y ,y s k y )- 
Under Hi, the number of counts fiij expected in a measured sky pixel is given by : 

$=b+8-fa (i) 

where b and s represent the (flat) background level and source strength respectively, and the 
number of source counts in pixel from a source at (x s ky,y s ky) (see Boella et al. 1997a). The 
measured number of counts in pixel is Nij. 

By minimizing the natural logarithm of the likelihood C under Hi, given by ln(£-^ 1 ) = 
Si (j^ij^-if^ij 1 ) ~ ^IJ 1 )' w ith respect to its free scale parameters, b and s, one can derive the 
flux and flux uncertainty from s and its error for a putative source at position (x s k y , y s ky)- From op- 
timizations under Hi and Ho one can determine the likelihood ratio A defined as — 21n(£-^ /£-^ 1 ). 
This quantity is distributed as x 2 f° r 1 degree of freedom for a known source position (e.g., Eadie 
et al. 1971) and yields the source detection significance. The maximum-likelihood ratios evaluated 
at the grid points (x s k y ,y s ky) constitute the MLR map. 

The MLR map for the MECS TOOl imaging data (8" grid point spacing; no energy selection 
applied) is given in figure 2. The maximum value corresponds to a ~ 12<r detection significance, 
and is located within the cross section of the WFC and IPN error regions. From the map it is clear 
that multiple sources contribute to the observed excess. For instance, the ~ 12<r excess is elongated 
and suggestive of two or more close point sources. In order to determine the minimum number of 
point sources on a flat background that are required to satisfactorily describe the measured data, we 
applied a maximum-likelihood source-fitting procedure. In this procedure the data are described in 
terms of a flat background model and a number of sources each with a free position. One starts with 
a model composed of the flat background model and a single source at an arbitrary free position and 
optimizes the likelihood with respect to the background level, source strength and source-position 
parameters simultaneously (4 free parameters). Next one determines the likelihood improvement 
with respect to the zero hypothesis Ho that the data are described in terms of a flat background 
model only (1 free parameter). If the improvement exceeds the 3<r threshold for 3 (= 4 — 1) d.o.f. 
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one continues and compares the optimized likelihood of a model composed of a flat background and 
2 sources at arbitrary positions with the optimized likelihood for the previous model with 1 source 
at an arbitrary position. This process is repeated until the likelihood improvement is below the 
3d threshold for 3 d.o.f. At each iteration one obtains the optimum background scaling parameter 
and for each source its strength and position together with uncertainty estimates on each free 
parameter. 
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Fig. 1. — Time history of the burst itself as seen with WFC and GRBM, at a time resolution of 
1 s. The vertical bars near the left edge of each panel indicate the la errors in the rates. 
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Fig. 2. — Maximum-likelihood contour map of the MECS data from TOOL The contour levels 
start at the equivalent of a 3a detection and increment with steps of la (1 d.o.f.). The highest 
contour has a level of 12a. The X and Y coordinates are with respect to the instrument coordinate 
frame. The symbols refer to the fitted positions of the four significant point sources detected 
in the combined TOOl and T002 data. The upward pointed triangle indicates the afterglow of 
GRB 991014 (SAX J0651.0+1136). The light-shaded circle is the WFC error region of GRB 991014 
(99% confidence) and the dark-shaded area between the two parallel lines bound the IPN annulus 
(3a confidence) 



